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Computing determinants is invariably in demand from all sorts of mathematical domains and needs 
neither cynical advocacy nor does it lack motivation. In this semi-expository article we wish to illustrate 
a few techniques that should receive ample attention due to their simplicity, elegance and effective 
nature. 

Glossary of notations and conventions 

Throughout, almost all determinants have indices ranging in {0, 1, . . . ,n — 1}. Where there is no 
confusion, we omit to write them. 

mil = 1!2! • ■ ■ ml and (x)i = x(x + 1) • ■ • (x + i — 1) is the Pochhammer symbol. 
(I) = if either 6 < or b > a. 

Cti is the constant-term operator that extracts the coefficient of x°. Also Ct = Yi^o 
S n is the symmetric group of permutations on n letters. 
x„ = (x ,xi, . . . ,x„_i) and X„ = x xi ■ ■ -x n -i. 
V(x n ) is the Vandermionde determinant rij^OEj ~~ x i)- 
M 9 = ^£and [n]!, = [l] 9 [2] g .-.[n] 9 

The Gaussian polynomials (£) = ^, ^^ fc j, , for < k < n; = for k < or k > n. 

1 The Method of Condensation 

The idea originated with Charles L. Dodgson and the iterative form is propelled by Zeilberger [7], Here 
are some instances of its application. 

1.1 Example. Consider the n x n matrix of entries Let us compute the determinant. The first 

step is to generalize it (introduce two additional parameters) as Z n (a, b) — det ( l ^_ b a ) ■ The next step 
is an automated guess for its determinantal value 

2(2) TT ^ + 2 ^ ! = 2 ®T\( 2i + 2a )( i + b 
1 ' W (i + b)\(2t + 2a-by. 11 V b J\ b 

Proof. According to [7], the following is always true for any determinant (when i — > i + a: j — > j + b) 

Z n -i(a, b)Z n -i(a + 1,6 + 1) - Z n -\{a + 1, b)Z n -i(a, 6+1) 



Z n (a, b) = 



Z n _ 2 (a + 1,6+ 1) 



It remains to verify that (1.1.1) also satisfies this non-linear recurrence. Symbolic softwares are naturally 
adept at this task. The proof will then be complete once two initial cases are checked, say for n = 1 
and n = 2. □ 
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1.2 Example. In particular, setting a = b = implies that dct (^) = 2W. The Online Encyclopedia 

of Integer Sequences (OEIS) associates 2( 2 ) with the number of n x n binary matrices with no row sum 
greater than n — 1, or the number of labeled n-colorable graphs on n nodes. 

We introduce a generalization of formula (1.1.1) that is provable by Condensation. 



(1.2.1) 



det 



ri + x 

i + v 



-(3) Y[ 



ri + x\ I 1 + y y 



i=0 



y 



v 



Once again, the specialization x = y = yields a pretty identity: det (") = A 2 \ Equation (1.2.1) 
also shows something interesting: the following determinant is independent of x. 



dct 



qi + x 
j . 



1.3 Example. The number of plane partitions that fit in an a x b x c box is enumerated by 



D n (a, b) := dct 



i + j + a + b 
i + a 



c-l 



z-l 6-1 0-1 . 



=nnn 

■;3=0 k=0 i=0 i=0 



i+j+k+2 
i+j + k + 1' 



Again, this formula is provable by the Condensation method. This determinantal representation is 
recorded here due to the apparent symmetry of the matrix on the left-hand side. 



2 Buy one, get three free! 

In the last example, set a = b = 0, and write D n instead for D„(0,0). We will calculate D n in three 
different ways, each leading to a different interpretation. 

2.1 A la Andrews 

The present method has been utilized by George Andrews [2] in his proof of the Mills- Robbins-Rumsey 
determinant [6] 

(2.1.1) A n :=dct 

The main essence was to exhibit a triangular matrix T n so that A n T n is triangular too. In the same 
spirit, let us implement the Vandermonde-Chu identity 

( V) = ^Q{i-k)' 

which follows from an elementary combinatorial argument. Therefore 




since in the latter the matrices are triangular with l's on the main diagonal. 



V+j + i 
2j-i 



n 



r (A^ + 2^ + 2) ^ (i A^ + 2^ + §) 
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2.2 Constant-term identities 

This method goes back, to at least, Zeilberger's paper [8]. For the current application, begin as follows: 

=Ct i (xi(l + xr 1 y+j)=Ct i (l + x i ) i (l + xr 1 y. 
Therefore 

n-l 

D n = Ct J] (1 + Xi y det [(1 + x-y] 

i=0 
n-l 

= CtJl{l + x i ) i '[l(xJ 1 -xr 1 ) 

= ±ci 1 [[(x j -x i )l[{xj 1 -xT 1 ) 

j>i j>i 
lli=0 x » j>j 

Since IIj>i( x j — ^i) i s ^ ne Vandermonde determinant V(x), we gather the coefficient of X™ -1 in V 2 (x n ) 
to be (— l)( 2 )n!. This, of course, reproves a trivial special case of Dyson's conjecture that the constant 
term in Ylj^ (l - fj-) equals n\. 

2.2.1 Remark. The even powers (in particular, the square) of the Vandermonde determinant play a 
crucial role in the quantum Hall effect phenomena via Laughlin's wave function ansatz. Determining 
the coefficients of their expansion, as symmetric functions, in terms of Schur functions has triggered 
considerable interest among physicists (see [3] and references therein). At the moment of this writing, 
extracting these numbers in general is an open problem. 

2.3 Multilinearity of the determinant 

Define the determinant B n = det[(i + For the third computation, watch an online video of a 
lively lecture delivered by Zcilberger [10]. For the purpose at hand, let us make use of Euler's formula 
to! = / °° t m e~ l dt. Thus 



B n = det 



/>oo 

/ x 1 ^ e~ Xi dxj 
Jo 



= L y e '^ X3V ^Jl x jj dXn 

= ^E / (e-^Kx-irnV'W 

' tt£S„ JR + \ 1=0 / 

= - f e-2>'V 2 (x„)dX n . 

Observe that B n = det !(*+■?)] = (n- 1)!! 2 £>„ = (n-l)!! 2 . Hence, we obtain the Selberg-type integral 

f e -S^y 2 (x„)dX„ = n!!(n-l)!!. 
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2.4 Remark. In summary, the above ideas combined with Dodgson's Condensation suggest simple 
techniques that are efficient in proving more general evaluations, such as 



det 



% + j + a 

i + (3 



c*n(i 



/ XZe-^^V 2 ^(^ n )dX n . 



3 More Examples 

The Delannoy numbers D(i,j) count walks from (0, 0) to composed of unit steps East, North, and 
North-East. These are given by J2 k >o (D GO 2 *- 

3.1 Example. The below determinant was conjectured by Bacher and Krattenthaler [5] provided a 
proof. We offer a simpler argument using Constant-term identities. For alternative automatic proof see 
Koutschan [4]. 

min(2i,2j) 



K n := det 



E 

k=0 



4(3) ft ™ 
AA (2i)!3 ■ 



Proof. Start with J2 k >o (fcHfc) 2 * = Ct i{ 1 + 2^)^(1 + xT x ) 2 K Then 



K n = CtH(l + 2x,f l det [(1 + x- 1 ) 



i=0 



= -^Ct JJfo + - a* - z 2 ) det [(1 + xr 1 ) 2 *] 

j>i 
n-1 

= S)ctl[(x i +x 2 ydet [(l + xr 1 )*] 
= 4© det [ctfo + x 2 ) l (l + x^fi 



= 4(3) det 
= 4(2) det 



ct i (^(i + a;i ) 2j+i ) 

2j + i 
2j-i 



Formula (1) in [1], with x = y, a = 2, computes the last determinant to produce the desired value. □ 

3.2 Example. Some play time with variants of Example 3.1 lead us to conjecture the determinant 
L n given below. The authors are grateful to Tiago Dinnis Da Fonseca (Univ. Montreal) who allowed 
his proof to be included in this article. It is rather amusing how L n turns out to be equivalent to 
(2.1.1). By association, its evaluation is as simple or as complicated. For alternative automatic proof 
see Koutschan [4]. 



L n := det 



min(2i,2j) 

E 



= 1fi(S )T T 1 M!M!(!_±_) 

Al (4»)!2(4i + l) 



»=o 
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Proof. It is evident that J2 k >o Ck) (ft) 4 * = ^( X + 2x i? l ( l + ^ 1 ? 3 ■ Th en 

n-l 

L n = Ct Yl (1 + 2x,) 2j det [(1 + 2xr l fi] 

i=0 

= ^-CiY[(xj + x) - a* - xj) det [(1 + 2xr 1 fi] 

j>i 
n-l 

= 4$)ct n (Xi + AY det [(1 + air 1 ) 2 '] 

n-l 

= &)& n + xir^) n^ 1 + ^ 2 - ^ x 

i=0 j>i 

= 16(2) Ct [] 0* + x iY det [fc rl + x i 2 Y] 

i=0 

)(2) det [ctfo + a; 2 ) 4 ^ 1 + , T r 2 )J 



16V 



= 16(2) det [Ct, (^"^'(l + Xi) i+j )] 
= 16(2) det 



j + i 
2j-i 



Alas! this is the familiar Mills- Robbins-Rumsey determinant (2.1.1), specialized at /U = 0, 



det 



j + i 
2j - i 



n-l 



(2z + 2),(2i + f)i_i 



2 "S Wi^+Di-i 



Using (y)j = ^rr^, the right-hand side routinely transforms to the required product expression. □ 
4 The Holonomic Ansatz 

Once more, this automatic method for the evaluation of determinants is a creation of Zeilberger [9]. 
We are grateful to Christoph Koutschan (RISC-Linz) for allowing his proofs, which exploit the present 
method of holonomic ansatz, to be made available [4]. So far, this is the only technique providing a 
proof to the following ex-conjecture. 



4.1 Example. We have the determinantal evaluation 



det 



min(3i.3j) 

E 



3j 



3 k 



9 8 (2) TT (^MAMlMll! 



6 2 (2) n n 



i=l j=l 



(12j - 5)(12j - ll)(4j + l)(4j - 1) 
(6j + l)(6j-5)(3j-l)(3j-l) ■ 
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5 Quantum Analogues 

After the initial release of the present note and its conjectures, the authors received instant feedbacks. 
They arc grateful to Johann Cigler (Univ. Vienna) for the permission to incorporate his conjectures 
and comments on certain g-analogues. Although the proofs might be carried out by employing the 
methods from the earlier sections, we are not pursuing the details. The reader is invited to do so. 



5.1 Example. A ^-analogue of formula (1.2.1) is given by 



det 



ri + x 

j + y 



i(y-x) 



, , n-1 U'll r ("+ x ) 

[r]Fn« (r - i)(S) - q 



i=0 



m + y v ) c 



5.2 Example. The determinants 

min(i,j) 



det 



E 



fc=0 



kj q \k 



=0 \ / q \ / q 



and 



det 



min(ij) , , N / , N fe 

E 



n( i+ ^) 



n-1 



i=0 



both follow from the general formula 
(5.2.1) det 



min(i,j) 

^ a(i,k)b(j,k)f(k) 

k=0 



n-1 



]J a(i,i)b(i,i)f(i). 



i=0 



The identity (5.2.1) results from an immediate application of the method outlined in Section 2.1. 
5.3 Example. The determinant 

min(ri,j) 



det 



n-1 



j=l »=0 



follows from 



det 



min(ri,j) 

E 



/•z 



= det 



min(ri,j) 

E 



rz 



n-1 



Therefore, it suffices to prove that 



det 



min(ri ,j) 



ri 



n—3 



i=0 
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